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We study the kinematics of nonlinear resonance broadening of interacting Rossby waves as mod- 
elled by the Charney-Hasegawa-Mima equation on a biperiodic domain. We focus on the set of 
wave modes which can interact quasi-resonantly at a particular level of resonance broadening and 
aim to characterise how the structure of this set changes as the level of resonance broadening is 
varied. The commonly held view that resonance broadening can be thought of as a thickening of the 
resonant manifold is misleading. We show that in fact the set of modes corresponding to a single 
quasi-resonant triad has a nontrivial structure and that its area in fact diverges for a finite degree 
of broadening. We also study the connectivity of the network of modes which is generated when 
quasi-resonant triads share common modes. This network has been argued to form the backbone 
for energy transfer in Rossby wave turbulence. We show that this network undergoes a percolation 
transition when the level of resonance broadening exceeds a critical value. Below this critical value, 
the largest connected component of the quasi-resonant network contains a negligible fraction of the 
total number of modes in the system whereas above this critical value a finite fraction of the total 
number of modes in the system are contained in the largest connected component. We argue that 
this percolation transition should correspond to the transition to turbulence in the system. 



I. INTRODUCTION 

The phenomenon of dispersive wave propagation is 
fundamental to our understanding of a wide variety of 
spatially extended physical systems. In such systems, 
the frequency, cjk, of each wave mode is a nonlinear func- 
tion of its wave- vector, k [1]. Examples include gravity- 
capillary waves on fluid interfaces [2], fiexural waves in 
thin elastic plates [3] , drift waves in strongly magnetized 
plasmas [4] and Rossby waves in planetary oceans and 
atmospheres [5]. In this article, we will be interested 
in Rossby waves as modeled by the Charney-Hasegawa- 
Mima (CHM) equation [4] on the ^^-plane: 



dt {A^ - F^) + Pd:,^ + J A^] = 0. 



(1) 



This is the simplest two-dimensional model of the large 
scale dynamics of a shallow layer of fluid on the surface 
of a strongly rotating spehere. The surface of the sphere 
is approximated locally by a plane, x G M^, with x vary- 
ing in the longitudinal (meridional) direction and the y 
varying in the latitudinal (zonal) direction. The field 
'0(x, t) is the geopotential height, f3 is the Coriolis pa- 
rameter measuring the variation of the Coriolis force with 
latitude, F is the inverse of the square of the deforma- 
tion radius and J [/, g] denotes the Jacobian of two func- 
tions, / and g. This equation admits harmonic solutions, 
V^(x, t) = Re [Ak exp(zk • x - iui^t)] with k G M^. These 
solutions, known as Rossby waves, have the anisotropic 
dispersion relation 



(2) 



Since Eq.(l) is nonlinear, modes with different wave- 
vectors couple together and exchange energy. If the non- 



linearity is weak, one finds that this energy exchange is 
generally quite slow and occurs most efficiently between 
groups of modes which are in resonance. For the CHM 
equation, such resonances involve three modes since the 
nonlinearity is quadratic. Four modes would be involved 
in the case of systems with cubic nonlinearity. Three 
wave- vectors (ki,k2,k3) satisfying the resonance condi- 
tions, 



ka ki +k2, 

cj(k3) - cj(ki) - cj(k2) = 0, 



(3) 



are referred to as a resonant triad. If one projects the 
spectral representation of the wave equation onto a reso- 
nant triad, one obtains a set of ordinary differential equa- 
tions for the coupled time evolution of the amplitudes of 
the constituent modes. Such systems of equations ap- 
peared as basic models of nonlinear mode coupling in a 
variety of physical systems including plasma physics [6], 
nonlinear optics [7] and oceanic internal waves [8]. An 
advantage of such models is that the equations of motion 
for a resonant triad are simple enough that explicit for- 
mulae can be obtained for both amplitudes and phases 
of the resonant modes [9-13]. 

A disadvantage, however, is that such triads are gen- 
erally not closed. Even if energy is initially mostly re- 
stricted to a single triad, other resonant modes can be 
generated which are not in the original triad. This pro- 
cess can repeat in a cascade-like fashion and result in the 
excitation of a large number of modes. If a large number 
of degrees of freedom are excited, a statistical description 
of energy transfer between modes is preferable. Such a 
description is provided by the theory of wave turbulence 
[2, 14]. This theory provides a kinetic description of en- 
ergy transfer in ensembles of weakly interacting disper- 
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sive waves in which conserved quantities are redistributed 
along the resonant manifolds. See [15, 16] for a review. 

For an infinite system, in which wave modes are in- 
dexed by a continuous wave- vector, the theory of weakly 
nonlinear wave turbulence becomes asymptotically exact 
in the weakly nonlinear limit. For finite sized systems, 
in which the wave modes are indexed by a discrete wave- 
vector, some subtleties arise. The simplest case, which is 
particularly relevant to numerical studies of wave turbu- 
lence, is a bi-periodic box. In this case, k, is restricted to 
a periodic lattice with a minimum spacing, A A:, between 
modes. Modulo this spacing, the components of k must 
be integer valued. This is an issue because if the com- 
ponents of k are integers, then the resonance conditions, 
Eq.(3), become a problem of Diophantine analysis. Such 
problems typically have far fewer solutions than their 
real- valued counterparts and it is generally quite difficult 
to find them. A complete enumeration of all solutions 
for the case of Eq.(2) with F = was recently provided 
in [17]. This sparseness of solutions means that, in dis- 
crete systems, resonant triads can exist in isolation or in 
finite groups of triads known as resonant clusters. Two 
triads belong to the same cluster if they share at least 
one mode. The dynamics of small clusters consisting of 
two triads has been studied in considerable detail in [18]. 
Small clusters have attracted some interest in the con- 
text of atmospheric dynamics as a possible explanation 
of the unusual periods of certain observed atmospheric 
oscillations [19]. Depending on the dispersion relation, 
there may or may not exist large clusters capable of dis- 
tributing energy over a large range of scales in a discrete 
system. In the case of Rossby waves on a sphere with 
infinite deformation radius, such a large cluster does ex- 
ist [20]. On the other hand, for capillary waves there 
are no exactly resonant triads at all [21]. For the dis- 
persion relation Eq.(2), numerical explorations indicate 
that for general values of F large exactly resonant clusters 
are rare. Thus, in discrete systems, it is often necessary 
to rely on approximate resonances to account for energy 
transfer. 

Approximate resonance is possible due to the phe- 
nomenon known as nonlinear resonance broadening. This 
is an effect whereby the frequency of a wave acquires a 
correction to its linear value which depends on the am- 
plitude (see Chaps. 14 & 15 of [1]). Triads which are 
not exactly in resonance can then interact at finite am- 
plitude if the frequency mismatch is less than this cor- 
rection. Such triads are known as quasi- resonant triads 
and satisfy the broadened resonance conditions 

rk3 = k, + k„ 

1 \uj{ks) - Uj{k^) - Uj{k2)\ < S ^ ^ 

where ^ is a characteristic value for the resonance broad- 
ening taken to be positive. Although Eq.(4) provides 
only a kinematic description of resonance broadening, 
the analogous dynamical effect can be very strikingly vi- 
sualised in linear stability analyses of weakly nonlinear 
waves [22, 23] where it is found that the set of unstable 



perturbations lie in a neighbourhood around the set of ex- 
actly resonant perturbations. This set of quasi-resonant 
modes is often pictured as a "thickened" or broadened 
version of the exactly resonant manifold. For weakly non- 
linear systems, this broadening is expected to be small 
since amplitudes are small. It may nevertheless be large 
enough to overcome frequency mismatches which arise 
when wave-vectors are restricted to a discrete grid and 
prevent discreteness from impeding the cascade of en- 
ergy. A striking example of this effect is observed for 
capillary wave turbulence in a biperiodic box. For this 
system, since there are no exact resonances, as the nonlin- 
earity is decreased the resonance broadening eventually 
becomes smaller than the frequency mismatches due to 
the grid. Direct numerical simulations illustrate that the 
cascade of energy to small scales stops entirely when the 
level of nonlinearity gets sufficiently small leading to the 
phenomenon of "frozen turbulence" [24-26]. 

The interplay between exactly resonant and quasi- 
resonant clusters means that wave turbulence in discrete 
systems is nowadays believed to exhibit several regimes. 
If the typical resonance broadening, ^, is small enough 
that effectively only exactly resonant clusters can inter- 
act, the dynamics are referred to as discrete wave turbu- 
lence [20, 27, 28]. If S is larger than the typical spacing 
between modes then effectively all triads can interact at 
least quasi-resonantly and the classical statistical theory 
is expected to be valid. In between is a regime consisting 
of a mixture of exactly resonant and quasi-resonant clus- 
ters which has been termed mesoscopic wave turbulence 
[28, 29]. In this intermediate regime, it has been sug- 
gested [30] that forced systems could exhibit some aspects 
of self-organised criticality. This suggestion is motivated 
by the idea that the forcing will cause the characteristic 
value of S to increase until it is large enough for a large 
quasi-resonant cluster to form which will then facilitate 
an "avalanche" of energy transfer to the dissipation scale 
thereby reducing wave amplitudes and the corresponding 
value of S. 

In this paper we develop a kinematic concept of criti- 
cality in quasi- resonant interactions in the Eq.(l). Specif- 
ically, we address the question of how a large quasi- 
resonant cluster emerges in the CHM equation as 5 is 
increased. Inspired by the theory of percolation on ran- 
dom networks [31], we take "large cluster" to mean a 
cluster that consists of a finite fraction of all modes in 
the system. We begin by analytically characterising the 
shape of the quasi-resonant set defined by Eq.(4) as a 
function of S for a single triad in Sec. II. By expressing 
the boundary of the quasi-resonant set in terms of the 
intersection of a pair of quadratic forms, we find some 
surprises. In particular, we find that the area of the set 
diverges at a finite value of S illustrating that the common 
perception of the quasi-resonant set as a "thickened" ver- 
sion of the exact resonant manifold is potentially quite 
misleading. In Sec. Ill we numerically construct the 
set of quasi-resonant clusters as a function of S for vari- 
ous system sizes. We show that a percolation transition 
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occurs at a critical value, of the resonance broaden- 
ing as S is increased. At this critical value, the size of 
the largest cluster rapidly goes from containing a negli- 
gible fraction of the modes in the system to containing 
a finite fraction of them. The value of decreases as 
the inverse cube of the system size, a fact which we trace 
to quasi- resonant interactions between small scale merid- 
ional modes and large scale zonal modes. We finish with 
a short summary and discussion about what conclusions 
can be drawn about Rossby wave turbulence from our 
results. 



II. CHARACTERISATION OF THE 
QUASI-RESONANT SET FOR A SINGLE TRIAD 

In what follows, we shall take ks to be fixed with k2 = 
ks — ki. The ^-detuned quasi-resonant set of k3 is the 
set of modes, ki, which satisfy the inequality 



|a;(k3)-a;(ki)-a;(k3-ki)| 



(5) 



This section is devoted to determining the structure of 
this set as a function of the detuning, S. The boundaries 
of this set are given by the pair of curves 

uj(ks)-uj(ki)-uj(ks-ki) = S, (6) 
uj(ks) - uj(ki) - coiks - ki) = -S. (7) 

We begin by finding these curves. Clearly it suffices to 
solve Eq.(6) since the second boundary can be obtained 
from this by setting S ^ —S. To fix notation, let us write 

ks = (p, q) , 
ki = (r, s) , 

7 2 2,2 

k = p -\- q . 

For the CHM dispersion relation, Eq.(2), the boundary 
of the quasi-resonant set, Eq.(6), then corresponds to the 
curve in the (x, y) plane implicitly defined by 



(3r 



f3{p-r) 



J^2^p^ ^2^g2^JP ' _ ^)2 + _ ^ ^ 



s. 



Subsequent formulae will be more compact if we shift the 
origin to the centre of symmetry of the curve, by defining 
X = r — p/2, y = s — q/2. Also, we rescale 6 by setting 
P = I from here on. The curve we wish to study is 
therefore 



P 



k^^F ' (x+p/2)2 + (i/ + g/2)2 + F 

X — p/2 

~ {x-p/2f + {y-q/2f+F 



(8) 



Gathering these terms together with a common denomi- 
nator we obtain the quartic curve c(x, ^) = 0, where 




FIG. 1. Solutions of Eq.(6) for different values of 5 with no 
self-intersection. Here ks = (cos ^, sin ^) with = 7r/6 and 
F = 1/5. The values of ^ -0.721688 and 62 ^ 1.20281 are 
obtained from Eqs.(ll) and (13) respectively. 



and the coefficients are given by 

ai = 1 + 4F) [3k^p - + F) + 4F) 5] , 
a2 = -p- + F) S, 

as = -I W + + 6F) - + - ^2 _ ^] ^ 

a4=l [P{P^ + 3*/' - 2F) - {e + F){p' - q2 + 4F) 5] , 
as = 2q[q^ + F- p{k^ + F) 5] . 

If we introduce new variables u = x'^,v = y'^ and w = xy, 
then it becomes clear that the boundary curve, Eq.(8), 
corresponds to the intersection of 2 quadratic surfaces 



a2{u + v)"^ + azu + a^v + ai = a^w , 



(10) 



W 



uv . 



c(x, y) = ai-\- a2{x^ + ^^)^ + a^x^ + a4y^ - a^xy, (9) 



Notice that the surface w'^ = uv is di cone, which is sin- 
gular at the origin x = y = 0. We have performed a full 
analysis of the properties of the intersection curves for 
general values of the parameters p, g, F and S. This is a 
technical exercise which is not very illuminating. In the 
interests of clarity, we will restrict ourselves here to illus- 
trating the essential qualitative features of these curves 
accompanied by some explicit examples. 



A. The curve typically exists only for a finite range 

of S, 



4 




FIG. 2. Solutions of Eq.(6) for different values of S which self- 
intersect Sit 5 = 82. Here ks = (cos^, sin^) with = Stt/T and 
F = 1/5. The values of ^1 ^ -0.185434 and 62 ^ 0.309057 
are obtained from Eq.(ll) and Eq.(13) respectively. 



This is evident from studying Eq.(8) as ^ ^ ±00. The 
LHS can only diverge if one of the denominators diverges. 
This is impossible if F > 0. Hence if F > 0, solutions 
to Eq.(8) must cease to exist if S gets too large by ab- 
solute value. In the exceptional case F = 0, solutions 
indeed exist for all values of 5 and become localised in 
the neighbourhood of =b(p/2,g/2) a,s 6 ^ ±00. 



The curve is bounded except at single critical 
value of 6. 



Examining Eq.(lO), it is clear that since u and v are 
both positive, it is generally not possible to balance the 
quadratic terms if + f ^ oc. Hence the curve, if 
it exists is bounded. The single exception to this occurs 
when the coefficient of the quadratic terms vanishes. The 
curve may therefore diverge at the single critical value of 
detuning given by: 



P 



(11) 



Note that this corresponds to 6 = uj{hs)- At ^ = ^1, some 
algebra shows that the curve is given by the hyperbola 



X -\-2-xy -y 
P 



F). 



(12) 



It follows that the area of the quasi- resonant set diverges 
as 6 tends to ^1 from below. This result illustrates that 
the common conception that the quasi-resonant set looks 



like a "thickened" version of the exact resonant manifold 
is a misconception. The special case p ^ corresponding 
to the case of ka becoming zonal is discussed separately 
below. 



The curve may self- intersect only at a single 
critical value 6. 



From Eq.(lO), self-intersection is only possible if the 
curve passes through (0, 0) in the v) plane. This can 
only happen if the coefficient ai = 0. Thus we identify 
a second critical value of S where self-intersection may 
occur: 



3k^p 



(13) 



Note that ai = does not necessarily mean the curve self- 
intersects. It is possible that at J = ^2, the curve reduces 
to a single point. To probe what happens at (5 = (^2, we 
consider the surface z = c{x^y) defined by Eq.(9) when 
S = S2 m the neighbourhood of the origin. Calculation 
of the partial derivatives indicate that this surface has a 
critical point at the origin. After some tedious algebra, 
we find that the determinant of the matrix of second 
derivatives is 

4(P + F)\k^ + 4(g^ - 3p^)F) 

If A > then the critical point at (0, 0) is a maximum 
or a minimum. The curve c(x, y) = is then an isolated 
point. On the other hand, if A < the critical point is 
a saddle and the curve c(x, y) = has a self-intersection 
at (0,0). The condition for self-intersection is therefore 



A:^-h4(g2 -3/)F > 0. 



(14) 



We note that for F = we always have a self-intersection. 

These qualitative features of the boundary of the quasi- 
resonant set are illustrated graphically in Figs. 1 and 2 
which show the shape of the curve for different values 
of S. Generic parameters, specified in the captions, have 
been chosen with no particular symmetries. Hence the 
curves shown in these figures are representative of what 
is found from the complete analysis of Eq.(8). Fig. 1 
shows an example in which no self-intersection occurs at 
S = 62, while Fig. 2 shows an example in which a self- 
intersection occurs. The hyperbolic curves identified at 
S = Si are clearly visible in both cases. 

A couple of special cases are worth noting: 

(i) p = 1, g = and F = 

This corresponds to a meridional mode. In this 
case, as = so the intersection of quadratic forms 
given by Eq.(lO) lies entirely in the w = plane 
and reduces to 

-(1 + 5){u + vf - 1(1 - - ^;) + 1(3 - <5) = 0. (15) 
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2 1 1 2 2 1 1 2 2 1 1 2 



{d)S = 2 (e)5 = 3 (f)(5 = 4 

FIG. 3. Shaded regions correspond to the resonant set defined by Eq.(4) with ka = (1,0), F = and /3 = 1 for different values 
of S. The exactly resonant manifold is the red solid line. 




2 1 1 2 2 1 1 2 2 1 1 2 



(a)(5 = (b)(5 = l/10 (c)(5 = l/2 

FIG. 4. Shaded regions correspond to the resonant set defined by Eq.(4) with ks = (0,1), F = and /3 = 1 for different 
values of 6. The exactly resonant manifold is the red solid line and corresponds to the two axes in this case. As explained in 
the remark below equation (16), this exactly resonant manifold consists of non-generic triads, i.e., triads where one or more 
interaction coefficients are identically zero. The case F = is chosen for simplicity only: for F > 0, the resonant sets have the 
same qualitative shape as in the case F = 0. 



We can immediately identify the critical points. 
The curve self-intersects at (5 = 3 and diverges at 
5 = —1. The point J = 1 is also noteworthy, being 
the complementary boundary to the divergent case. 
For this value of 6 the curve is a perfect circle. 



(ii) p = 0, g = 1 and F = 

This corresponds to a zonal mode. In this case the 
curve simplifies to 

5 (u v)'^ ^ - S (u - v) ^ ^ S = -2w. (16) 
2 16 
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The only special value of 5 for this case is ^ = 0. We 
then recover the exact resonant manifold of a zonal 
mode, xy = 0. This consists of the two coordinate 
axes. It is now less surprising that the boundary of 
the quasi-resonant set can diverge for finite 6 once 
we appreciate that the exact resonant manifold of 
a zonal mode is unbounded. The divergence of the 
boundary of the quasi-resonant set of the non-zonal 
modes in some sense reflects the presence of this 
structure in the dispersion relation. 

Remark. Notice that the exactly resonant manifold in 
case (ii) consists of non-generic triads, i.e., triads where 
one or more interaction coefficients are identically zero. A 
point on the x-axis corresponds to a so-called catalytic in- 
teraction (|ki| = |k2|), for which the zonal mode ks does 
not change its energy but influences the energy exchange 
between the other two modes in the triad. A point on the 
y-3x.is corresponds to a 'spurious' triad, formed by purely 
zonal modes, which do not interact at all: the interac- 
tion coeflicients are all identically zero because the three 
modes are collinear. In our computations of the network 
of quasi-resonant modes (Section III), non-generic triads 
are discarded from the start. 

Having given a fairly complete qualitative description 
of the behaviour of the curves which deflne the bound- 
aries of the resonant sets as the detuning is varied, it 
now remains to assemble these boundary curves for pos- 
tive and negative values of 5 to determine the interior 
and exterior of the quasi-resonant set. This is again best 
accomplished by illustration. Figs. 3 and 4 illustrate 
the quasi-resonant set for the two special cases discussed 
above for a range of increasing values of S. The shaded 
areas in these flgures correspond to the quasi-resonant 
sets. The exact resonant curve is also shown for reference. 
Fig. 3 illustrates one of the key points of this article: the 
common picture of the quasi-resonant set as a thickened 
version of the exact resonant manifold is appropriate only 
for small values of the broadening (eg Fig. 3(b)). One 
might counter this with the observation that it is only 
in the weakly nonlinear regime that it makes sense to be 
discussing quasi-resonant interactions in the flrst place 
and in this regime, the broadening is necessarily small. 
Eq.(ll) tells us, however, that no matter how small the 
broadening, there are always modes close to the zonal 
axis, whose quasi-resonant set diverges. 




FIG. 5. The normalised total number of clusters, plotted 
as a function of S for different system sizes, M. Here, we have 
fixed F = /3 = 1. 



M=100 — K- 
M=200 — B- 
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M=1000 
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FIG. 6. The density of the largest cluster, p, as a function 
of resonance broadening 6 for different system sizes, M. Here 
we have fixed F = /3 = 1. The percolation threshold, 
decreases as a function of system size. The inset shows the 
same data plotted as a function of 6 . The solid line is a 
fit of Eq.(19) to the data. 



III. STRUCTURE OF THE NETWORK OF 
QUASI-RESONANT MODES 

In a turbulent system many modes are excited. The 
fact that a mode can be a member of more than one 
quasi-resonant set leads to overlap between sets. This 
allows modes to join together to form a network of quasi- 
resonant clusters analogous to the exactly resonant clus- 
ters which have been extensively studied in the literature 
to date. In this section we study the structure of this 



quasi-resonant network as the typical amount of broad- 
ening, (5, in the system is varied. We consider a finite sys- 
tem of size 27r x 27r and truncate the wavenumber space 
such that there are M modes in the and ky direc- 
tions, modes in all. The spacing between modes is 
2it/M . The quasi-resonant clusters were obtained by an 
exhaustive numerical search speeded up by incorporat- 
ing symmetries of the triads. As explained in the remark 
below equation (16), we have eliminated from our list of 
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triads the so-called non-generic triads, i.e., the triads for 
which one or more interaction coefficients are zero. This 
includes triads formed out of collinear modes and triads 
where any two wave- vectors have the same modulus. 

We ffist measure the total number of clusters nor- 
malised by the maximum number of possible isolated 
clusters which is approximately This is shown 

as a function of S for several different system sizes, M, 
in Fig. 5. For each system size, as S increases, the total 
number of clusters ffist increases, then reaches a maxi- 
mum at a particular value of (5, which we shall denote 
by S^^ and then decreases. The value, at which the 
maximum occurs decreases as the system size, M, is in- 
creased. In order to understand these results one must 
realise that for the CHM dispersion relation, the num- 
ber of exactly resonant clusters which live on the grid of 
discrete wavevectors is rather small with most discrete 
triads necessarily exhibiting a small amount of detun- 
ing enforced by the geometry. The initial growth in the 
number of clusters is explained by the addition of isolated 
triads to the list of quasi-resonant clusters as the broad- 
ening grows large enough to incorporate the geometrical 
detuning generated by the grid. While one would expect 
this rate of increase to slow down as clusters start join- 
ing to form larger clusters, the attainment of a maximum 
and subsequent decrease in the total number of clusters 
becomes much easier to understand when we measure the 
size of largest cluster, p, which is simply the number of 
modes in the largest cluster normalised by the system 
size, M^. This is shown in Fig. 6. We see that as S ap- 
proaches the largest cluster quickly makes a transition 
from including a negligible fraction of the total number 
of modes in the system to including almost all of them. 
Thus the network of quasi-resonant modes undergoes a 
percolation transition ai S = S^. While the results shown 
in Figs. 5 and 6 were obtained for /3 = F = 1, we checked 
that the results are not very sensitive to this choice. Al- 
though this is an entirely kinematic study, we expect that 
the emergence of a percolating quasi-resonant cluster is 
the explanation for the dynamical transition from dis- 
crete wave turbulence to Kolmogorov turbulence in the 
CHM equation as the nonlinearity of the system is in- 
creased. 

The decrease of the percolation threshold, S^^ as the 
system size, M, increases is clearly due to the fact that 
as M increases, the spacing between discrete modes de- 
creases. The amount of broadening required to overcome 
the geometrical detuning is therefore smaller and hence 
clusters can form more easily. We would like to quantify 
this decrease. We have already learned that the area of 
the quasiresonant set of any particular mode, k, diverges 
at a finite value of 5 = cj(k) as indicated in Eq.(ll). It 
seems plausible that as soon as the broadening is suffi- 
ciently large to allow any mode in the system to have 
a divergent quasi-resonant set, then the largest cluster 
must be of the size of the system. Thus one might esti- 




M 

FIG. 7. Percolation threshold, plotted against system size 
and for different values of F. Here, we have fixed /3 — 1. For 
large system sizes, exhibits the scaling ^ M~^'^ , where 
we measure ai = 3.00 (2.97,3.03). The brackets show the 
95% confidence interval. 



mate 

J* ^min|a;(k)|. (17) 

k 

By this argument, Eq.(2) tells us that for large systems 
we should see the scaling, 5^ ~ M~^. A numerical mea- 
surement of how (5* scales with system size for a range 
of values of F, is plotted in figure 7. We see that the 
value of F becomes redundant for large M as already 
mentioned. Each of the curves converge to the straight 
line indicated by the scaling 

(5* -M-^S (18) 

with exponent ai = 3.00 with a 95% confidence intervals 
if (2.97, 3.03), which was obtained by bias-corrected boot- 
strapping. The connected component therefore forms 
much more easily than the naive argument above would 
suggest. The origin of this scaling can be traced to 
the fact that zonal modes require very little detuning in 
order to interact with high k almost-meridional modes. 
Take F = for simplicity, although the following argu- 
ment works for any F > 0. Let us consider the largest- 
scale zonal mode in the system, k = (0, 1). We recall that 
the exactly resonant manifold xy = {) oi this zonal mode, 
figure 4(a), gives rise to non-generic triads which do not 
interact efficiently and therefore must be discarded. So 
we ask what is the minimum value of the detuning so 
that the quasiresonant set of the mode k = (0, 1) con- 
tains some new modes. From Fig. 4(b) it is clear that 
the first new mode to join the quasi-resonant set will be 
ki — (M, 0), corresponding io x = M^y = —1/2. By 
directly replacing these values of x^y into Eq.(16), we 
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FIG. 8. Color map in k-space showing the smallest value 
of broadening required for each mode to become a member 
of a cluster of any size. The darker regions are resilient to 
becoming members of clusters in the sense that large values 
of broadening are required. Lighter regions join clusters easily. 
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FIG. 9. Fitting the data in Fig. 6 to a standard phase tran- 
sition profile. The lefthand side of Eq.(19) is plotted as a 
function of S. The fit (solid line) is taken over the range 
[6^ : 2] with 6^ = 0.815. It has slope 2.03. which suggests a 
value of 2; ~ |. 
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16 



M. 



For large M, this has solution ~ M~^, in agreement 
with Fig. 7 and the inset of Fig. 6. This suggests that the 
percolation transition is driven by interactions between 
large scale zonal modes and small scale meridional modes. 
This is consistent with the known scale-nonlocality of 
wave turbulence in the CHM equation (see [32] and the 
references therein) and provides further evidence, albeit 
circumstantial, that the percolation transition is associ- 
ated with the onset of turbulence in the CHM model. 

In order to further illustrate this point, we can ask 
which modes have the greatest or least tendency to be- 
come part of a quasi-resonant cluster. This is done in 
figure (8), where we have coloured each mode according 
to the minimal amount of resonance broadening required 
for this mode to join a cluster of any size. Blue modes be- 
come active very easily whereas red modes are resilient 
to becoming part of any cluster. Rather than appear- 
ing homogeneous and random, we see the appearance of 
definite structure. We see a circular region containing 
modes with a strong propensity to join clusters including 
a narrow large-scale region of zonal modes accumulating 
near the kx = axis with very low interaction thresh- 
olds as expected from the discussion above. We also 
remark upon the group of large-scale meridional scales 
reluctant to form any quasi-resonant connections. The 
relative reluctance of large-scale meridional modes to ex- 
change energy has already been remarked up in the lit- 
erature and suggested as an explanation of the inherent 



anisotropy of Rossby waves turbulence. For example in 
[33] , a wave-turbulence boundary was computed by com- 
paring the CHM dispersion relation to the inverse of the 
eddy-turnover time. It is then argued that inside this 
region Rossby waves dominate, but with a frequency in- 
commensurate with that of the surrounding turbulence 
so that energy cannot penetrate into this region. The 
boundary thus obtained seems similar in position and 
shape to the dark region in Fig. 8. This anisotropic 
energy distribution at large scales has been documented 
in numerical simulations of Rossby wave turbulence (see 
[34] and the references therein). It is somewhat surpris- 
ing to see it emerging again here from purely kinematic 
considerations. 

Finally we might ask whether the density of the giant 
cluster illustrated in Fig. 6 shows the generic profile for 
a second order phase transition. 



piS) 







S<5, 



(19) 



and, if so, what is the exponent z. Notice that the system 
size M has been absorbed after appropriate rescalings of 
S and by the factor M^. Eq.(19) contains 3 adjustable 
parameters, (5*, c and z, which makes it difficult to un- 
ambiguously determine the exponent z. As pointed out 
in [35], if Eq.(19) holds, then 



d5 =z^^- 



(20) 



for S > S^. Plotting this quantity against S produces 
an easier fitting problem because the amplitude c, can- 
cels out, the fitting becomes linear and the value of 5^ 
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corresponds to the point at which the fitted hne crosses 
zero. The values of ^ would ideally be obtained inde- 
pendently from the values for p. In our case, this was not 
possible so we obtained them by locally interpolating the 
measured values of p using Mathematica and differenti- 
ating the result. The outcome of this analysis is shown 
in Fig. 9. We can see that a case can be made for a 
linear fit in the neighbourhood of S^. Taking = 0.815 
(the point at which the straight line fit obviously starts 
to fail), and fitting the data over the range [S^ : 2] gives 
the fit shown in the figure. The slope gives a value of 
z ~ ^. This is standard Landau value for the mean field 
theory of a second order phase transition of a scalar field. 
These results, while suggestive, are far from definitive. A 
more detailed numerical study in the vicinity of will be 
required before we can start putting estimates of uncer- 
tainty on these values. For the purposes of comparision, 
Eq.(19) with the best fit values of the parameters, is plot- 
ted with the original data in the inset of Fig. 6 (solid line). 



IV. CONCLUSIONS AND OUTLOOK 

To conclude, we have presented a kinematic analysis 
of the properties of quasi-resonant triads in the CHM 
equation. We described the analytic form of the quasi- 
resonant set defined by the quasi- resonance conditions, 
Eq. (4) as a function of the resonance broadening, S. We 
found that they have non-trivial geometric shape and are 
not well described as simply thickened versions of the 
exact resonant manifold as is commonly assumed. In 
particular, we found that the quasi- resonant set becomes 
unbounded for above a critical value of S and that this 
can occur for arbitrarily small values of S as we consider 



modes approaching the zonal axis. We then conducted 
an in-depth numerical study of the structure of quasi- 
resonant clusters as a function of 6 and identified a per- 
colation transition as 6 is increased. At the transition, 
a large cluster is formed which contains a finite fraction 
of all the modes in the system. For a system contain- 
ing modes, the value of the percolation threshold 
decreases as M~^. This scaling results from the ease 
with which large scale zonal modes interact with small 
scale meridional modes, a reflection of the nonlocality of 
Rossby wave turbulence. 

We speculate that the percolation transition corre- 
sponds dynamically to the transition from mesoscopic to 
classical wave turbulence. The fact that a percolation 
transition exists is consistent with earlier work on capil- 
lary waves [26]. In fact, we believe that this transition is 
not a special feature of the CHM dispersion relation and 
is generic [36]. Furthermore our results provide circum- 
stantial support for the sandpile picture of mesoscopic 
wave turbulence suggested in [30] since a small change in 
resonance broadening in the vicinity of can trigger a 
transition from a state which cannot support an energy 
cascade to one which can. In order to better understand 
these issues, we believe that it is important to move be- 
yond the kinematic picture of resonance broadening and 
attempt to devise methods of studying these effects dy- 
namically. 
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